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Wootters [PRL 80, 2245 (1998)] has derived a closed formula for the entanglement of formation 
(EOF) of an arbitrary mixed state in a system of two qubits. There is no known closed form 
expression for the EOF of an arbitrary mixed state in any system more complicated than two 
qubits. This paper, via a relatively straightforward generalization of Wootters' original derivation, 
obtains a closed form lower bound on the EOF of an arbitary mixed state of a system composed of 
a qubit and a qudit (a rf- level quantum system, with rf > 3). The derivation of the lower bound is 
detailed for a system composed of a qubit and a qutrit {d = 3); the generalization to (i > 3 then 
follows readily. 
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I. INTRODUCTION. 



As Wootters [1] discusses, the " entanglenient of formation" (defined below, hereinafter simply "EOF") is one of the 
more commonly employed measures of the entanglement of an arbitrary mixed state of a bipartite quantum system, 
i.e., a quantum system composed of two and only two quantum subsystems; for pure states of bipartite quantum 
I , systems, the von Neumann entropy (also defined below) is the standard measure of entanglement. Indeed Wootters 
J>. ' [2] , taking advantage of the fact that the von Neumann entropy of a pure state of two qubits is a function of a single 
, real parameter only, has derived a closed formula for the EOF of an arbitrary mixed state in a system of two qubits. 
T-H There is no known closed form expression for the EOF of an arbitrary mixed state in any bipartite system more 
complicated than two qubits. As elaborated below, however, the von Neumann entropy of a pure state composed of 
' ^ a qudit (any wavefunction is a linear combination of at most d orthonormal eigenfunctions) and a qubit (a qudit in 
the special case d — 2) also is a function of a single real parameter only. 

This paper therefore has sought to generalize Wootters' two qubit derivation [2], so as to obtain a closed formula 
for the EOF of an arbitrary mixed state of the qubit-qudit system. This attempt has not been successful, for reasons 
that will be manifest, but it has proved possible to derive a closed form lower bound on the EOF of an arbitrary 
qubit-qudit mixed state. Presenting this derivation is the primary objective of this paper. The derivation will be 
r-| ' detailed for an arbitrary mixed state of a qubit and a qutrit (a qudit in the special case d = 3); the generalization to 
the qubit-qudit (with d >3) case then follows readily. 
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II. THE QUBIT-QUTRIT SYSTEM. 



Let p denote an arbitrary density matrix of a bipartite quantum system S composed of a qubit and a qutrit. Unless 
otherwise stated, it is assumed throughout this paper that ^ p, i.e., that p does not fortuitously happen to represent 
a pure state. Then the mixed state density matrix p can be written in the form 

Af 

p=^Po^'o*i (1) 

a=l 

(i.e., can be "decomposed") in an infinite number of ways [3]. In Eq. (1): p is a 6x6 matrix; each p^t is a real number 
> 0; SqPq = 1; 2 < r < 6 is the rank of p, i.e., the number of non-zero eigenvalues of p; is > r and can be an 
arbitrarily large integer; the dagger denotes the adjoint; Dirac notation, which in this writer's opinion beclouds the 
matrix structure of the pertinent mathematics, deliberately has been avoided (as it will be during the remainder of 
this paper); and each is a normalized system S wave function, expressible as 

2 3 

i=i j=i 

where the Ui and Vj are orthonormal eigenfunctions of the qubit and qutrit respectively. In general the different ^q, 
forming any given decomposition (1) are not orthogonal to each other, but for each ^! a the coefficients a"- satisfy 
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J2K\' = H^c. = i. (3) 



According to the Schmidt decomposition theorem [4], for any Eq. (2) can be replaced by 

*„ = 5^c^«M, (4) 

k 

where: the are orthonormal hnear combinations of the Uj; the are orthonormal hnear combinations of the 
Vj; and the quantities c" are non- negative real numbers. As the notation indicates, all the u", V^, and are a- 
dependcnt. Since there arc at most two orthonormal for any a, there are at most two non-vanishing in the 
Schmidt expansion (4) for \1/q,, even though there can be three orthonormal u^. Correspondingly the normalization 
relation (3) reduces to 

E(c^)' = 1, (5) 
fe=i 

wherein the range of each is < < 1, of course. 

The von Neumann entropy associated with the qubit or qutrit subsystem reduced density matrices that can be 
constructed from the wave function of Eq. (2) is [1] 

^;(*„) = -;f^(c^flog2(c^)^ (6) 

fe=l 

where the are the coefficients in the expansion (4). Because of Eq. (5) -E(\l/a) actually is a function of a single 
parameter only, which conveniently can be chosen to be 

C«=2c?c^; (7) 

evidently < < 1. Eq. (5) then yields, for use in Eq. (6), 

(c?)' = ^(l + v/l^) (8) 



(C2f= 2(1-V^^) (9) 
where we now are choosing c" > C2 . Therewith we can rewrite Eq. (6) as 

= (10) 

with e{Ca) the expression obtained after substitution of Eqs. (8)-(9) into the right side of Eq. (6). 
The EOF of p is defined to be [1] 

N N 

E(p) = min J2 PaE{^a) = min ^ Pae{C„), (11) 

a=l a=l 

minimized over all possible decompositions of p of the form (1).^ Also define the similarly minimized quantity 

JV 

C(p) = min 

X]fa<^a- (12) 
Because [1] e{Ca) is a monotonically increasing convex function of in its range < Cq, < 1, 



^We shall not concern ourselves here with the difference between the minimization operation indicated in Eq. (11) and the 
infimum operation favored in Wootters' review [1]. 
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a a 



(13) 



implying [1] 

e[C{p)] < E(p). (14) 

In order to obtain a closed formula for E(/3) from Eq. (14), or even merely a bound on E(/9), it is necessary 
to establish an explicit connection between C„ and ^'q, well beyond the implicit connection inherent in Eq. (4). 
Wootters' ability to establish such a conncc;tion. via introduction of the "concurrence" [2], is a key feature of his 
successful derivation of a closed formula for the EOF of an arbitrary mixed state in the two qubit case. At this 
juncture, therefore, it is helpful to briefly review Wootters' definition and application of the concurrence. 



A. The Concurrence. Wootters' Two Qubit Derivation. 

Until further notice we are considering the two qubit system only, rather than the qubit-qutrit system S to which 
the preceding Eqs. (1)-(14) pertain, p now is a 4x4 matrix, not 6x6; the range of r in Eq. (1) now is 2 < r < 4; and 
the subscript j in Eqs. (2) and (3) now is permitted to range over the values 1,2 only, rather than the values 1,2,3. 
With these understandings Eqs. (1)-(14) are applicable to the two qubit system without further alteration. Then in 
effect the concurrence of any two qubit wave function is defined as [2] 

C7on(«'„)= |*ts4.*|, (15) 

wherein: is a column matrix whose components ^'ma, m = 1, 2, 3, 4, are respectively the coefficients a"]^, a"2, ^21' ^22 
appearing in the two qubit version of Eq. (2); the asterisk denotes the complex conjugate; the notation is consistent 
with Eq. (1), with 4"]^ the row matrix that is the transpose of ; and S is the matrix 



S = 



/ 0-1 
10 
10 

V -1 



(16) 



Evidently 

Con(*„) = |2(a?ia^2-«2i«?2)l- (17) 
The numerical coefficients ag- in Eq. (2) now form a square matrix 



to which the singular value decomposition theorem [5,6] 

A„ = [/„A„y„ (19) 
applies. In Eq. (19): Ua and Va are (a— dependent) unitary matrices; the diagonal matrix 

A.= (?2«); (20) 

and c",C2 are the coefficients appearing in Eq. (4) with (as previously specified) cf > C2- Eqs. (18)-(20) immediately 
imply the determinantal equality |det(Aa)| = |det(Aa)| , i.e.. 



111^22 "21"'12 I ~ "-l ''2 ■ 



(21) 



It follows that Con(4'a) is precisely equal to the parameter defined in Eq. (7). Thus the two qubit version of Eq. 
(12) becomes 



JV JV 



C{p) = mm J2 Pa \H^K\ = min^p„ *aS*a , (22) 



a=l a=l 
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where is the transpose of '^a- Evidently C(p) is the minimum average concurrence that any decomposition of p 

can attain. 

Starting from Eq. (22), Wootters [2] obtains his closed form expression for C(p) as follows (in effect). Eq. (1) 
shows the matrix p is not only Hermitian but also is positive semidefinite, meaning [10] that 0^ pQ > for any two 
qubit wave function 0. Thus p can be brought to diagonal form by a unitary transformation Up, 

4 

p = Up\Ul = ^ ^lsWswl (23) 

s=l 

wherein the /ig (/ii > /i2 > ^ ^ 0) ^re the eigenvalues of p and the Ws (components Wms, n^i = 1,2,3,4) 
comprising the columns of Up are a corresponding set of orthonormal eigenfunctions. Because p is a density matrix, 
i.e., has trace Tr(p) = 1, Eq. (23) has the form of Eq. (1), i.e., Eq. (23) provides a special decomposition of p. 
Introducing Wg = ,/JIl'Ws, and defining W to be the 4x4 square matrix whose elements are Wmsi Eq. (23) takes the 
form 

4 

p = ^WsWl = WW^ (24) 

wholly consistent with matrix notation. Similarly, introducing <1>q, = and defining * to be the 4 x rectangular 

matrix (4 rows, A/' columns) whose elements are ^ma, Eq. (1) becomes 

N 

p=Y,^cc^i = ^^\ (25) 

again wholly consistent with matrix notation. Furthermore, starting with any set of column matrices satisfying 



Eq. (25), defining = ^a/y (^a^a)) and remembering Tr(p) = 1, it can be seen that such a set satisfies Eq. 
(1). In other words not only does any decomposition of p provide a set satisfying Eq. (25), but also any set 
satisfying Eq. (25) provides a decomposition of p. 
Next let R be any 4xA^ matrix satisfying 

RR+ = I, (26) 

where I is the unit matrix. Eq. (26) implies that the 4 rows of R (but not necessarily its A'^ columns unless R happens 
to be square) form a set of iV-component orthonormal vectors. Then, as Wootters [2] observes, the matrix $ defined 

by 

* = WR (27) 

satisfies Eq. (25). Moreover it can be proved [3] that every $ satisfying Eq. (25) necessarily satisfies Eq. (27) for 
some R obeying Eq. (26). It follows that Eq. (22) can be replaced by 

N N 

C(/)) =min^ |(*S*)„a =min^|(RWSWR)ac« , (28) 

a— 1 OL—1 

minimized over all possible R satisfying Eq. (26). Because S defined by Eq. (18) is symmetric, the square matrix 
WSW appearing in Eq. (28) also is symmetric. Consequently, via an extension of the singular value decomposition 
theorem to symmetric matrices [9] , 

WSW =Uw^wUw, (29) 
where Uw is unitary and Aiy is a diagonal matrix whose diagonal elements 

Ai > A2 > A3 > A4 > (30) 
are the positive square roots of the assuredly non-negative real eigenvalues of the Hermitian matrix 

WSW(WSW)^ = WSWW^SW* = WSpSW*. (31) 
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It can be seen that the eigenvalues of WSpSW* coincide with the eigenvalues of W*WS/?S =p*SpS, as well as with 
the eigenvalues of pSp*S and /3^/^Sp*Sp^/^; the matrix p^/^Sp*Sp^/^, hke WSpSW*, is obviously Hermitian and 
positive semidefinitc. 

Substituting Eq. (29) into Eq. (28) fixes our attention on the matrix Q = Uw^, which like R is a 4 x iV rectangular 
matrix. Also since 



QQt = C/iyRRt[/t = UwUl = I 



(32) 



the four rows of Q, again like the four rows of R, form a set of N-component orthonormal vectors. Moreover to every 



Q satisfying Eq. (32) there corresponds an R = U^Q satisfying Eq. (26). Thus Eq. (28) simplifies to 



JV 



N 



C(p) = min^ (QAvyQ)ca =min^ 



a=l 



2 

ma 



m=l 



(33) 



minimized over all possible Q satisfying Eq. (32). In fact (as Wootters [2] observes), because of the absolute value 
signs the minimization in Eq. (33) can be restricted without loss of generality to matrices Q which not only satisfy 
Eq. (32) but for which also every Qi^ is real and > 0. 
The trace of the matrix QAvfQ obeys 



Tr(QAvKQ) 



N 4 



a=l m=l 



2 

ma 



N 



E^- Eq 



2 

ma 



m=l 



Ka=l 



N 

^E 

a=l 



2 

ma 



(34) 



Therefore, now restricting Q as described at the end of the preceding paragraph, Eqs. (33) and (34) imply 



mm 



N 



N 



N 



Al+A2^QL + A3EQ3a + A4EQ: 



< C(p), 



(35) 



a=l a— 1 a—1 

wherein the coefficients of A2, A3, A4 are complex numbers having absolute values < 1. It follows that if 

Ai - A2 - As - A4 > 0, (36) 

then the smallest possible value of the minimum on the left side of Eq. (35), which minimum furnishes a lower bound 
to C(p), is Ai — A2 — A3 — A4. In fact this minimum actually is as small as Ai — A2 — A3 — A4, because with the 
so-restricted 




(37) 



the sums J2a==i Qma ™ Eq. (35) attain the value -1 for m equal to each of 2, 3 and 4; with this same Q = Q>, 
moreover, the right side of Eq. (34), which by virtue of Eq. (33) furnishes an upper bound to C(p), also equals 
Ai - A2 - A3 - A4. Thus when Eq. (36) holds, Eqs. (33) - (35) imply 



C(p) = Ai - A2 - A3 - A4. 

On the other hand, when 

Ai - A2 - A3 - A4 < 0, 

it readily can be seen (again as Wootters [2] observes) that there always exist real angles 62,03,04 for which 



Ai + A2e"^''^ + Ase^*"^ + A4e 



0. 



Correspondingly when Eq. (39) holds 



C(p) = 



(38) 

(39) 
(40) 
(41) 



because, using Eq. (40) together with the matrix 
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Q-Q< = i 



1 


1 


1 


1 




















_gie4 


_gie4 


gie4 



(42) 



one sees that the right side of Eq. (34) and the left side of Eq. (35) equals zero. The matrix Q<, which also is 
restricted as described above, was given (in essence) by Wootters [2]. 
Wootters [2] has embodied Eqs. (38) and (41) in the single equation 

C(p) = max(0,Ai - A2 - A3 - A4). (43) 

Note, as Wootters [2] recognizes, that the matrix Q> is not the only 4x4 (restricted as described above) unitary 
Q which, after insertion into Eqs. (33) - (35), implies Eq. (38) when Eq. (36) holds. In particular if the 4x4 
matrix P is orthogonal, i.e., is real and satisfies PP = I, then Q = Q>P can serve, because (i) rr(PQ> Avi/Q>P) 
= rr(Q> AvyQ>PP) = rr(Q> Ah'Q>), and (ii) the expressions ArnQma t;he right side of Eq. (34) remain 

real for all a; it is additionally required only that each of these expressions be positive, as they are for Q = Q>, a 
requirement which therefore always is achievable for some range of P as P is continuously varied from the unit matrix 
I. There also is a range of so-restricted unitary 4x4 matrices, all of the form of Eq. (42), which imply Eq. (41) 
when Eq. (39) holds, because in this circumstance (assuming A4 ^ 0) there evidently is a range of angles 62,03,64 
consistent with Eq. (40). It is additionally evident that the matrices Q yielding Eq. (43) need not be square, i.e., 
that the number A'^of terms ^^a^a in a decomposition of p yielding the minimum possible average concurrence can be 
greater than 4; in particular, the real matrix P satisfying PP = I introduced earlier in this paragraph can be 4xA'', 
with N arbitrarily large. 

Actually the foregoing review of Wootters derivation of Eq. (43) has implicitly assumed that the rank r of p is 4. 
If any fXg in Eq. (23) happens to be zero, i.e., if r is less than 4, the corresponding column Ws of W is identically 
zero. But Eq. (27) makes the elements of $ completely independent of the rows of R corresponding to identically 
zero columns of W. It follows that when r equals 2 or 3 Eq. (27) can yield a 4" satisfying Eq. (25), i.e., can yield a 
decomposition of p, even though R does not fully satisfy Eq. (26). Indeed the only rows of R, as well as of Q = 
(it can be seen), which surely form a set of N-component orthonormal vectors when r equals 2 or 3 are those rows 
which do not correspond to the identically zero columns of W. It can be shown, however, as Wootters [2] remarks, 
that whenever a /i^ = the corresponding A,, = in Eq. (35); thus the possibly deviating elements of Q appear in 
Eq. (33 - (35) only as coefficients of these vanishing Ag, and cannot affect the validity of the inferences we have drawn 
from those equations. In particular these possible deviations of Q from strict compliance with Eq. (32) do not alter 
the conclusions that: (i) no decomposition of p can yield a smaller average concurrence ^^^PaCa than is given by 
Eq. (43), and (ii) there always is at least one decomposition which, via Eqs. (37) and (42) together with R = Uy^Q, 
actually does yield an average concurrence equal to the right side of Eq. (43). 



B. Derivation of the EOF Lower Bound 



In the qubit-qutrit system, to which we now return, the analog of the matrix on the right side of Eq. (18) for \E'a 
of Eq. (2) is 

( "11 "12 '^la ] (A A) 

\ "21 "22 "23 / 

which no longer is square. The singular value decomposition theorem is applicable to non-square matrices [5,6], but 
for the purpose of computing the determinants of matrices, as was done in deriving Eq. (21), it is necessary [7] to 
work with square matrices. Therefore the qubit-qutrit analog of Eq. (18) will be taken to be 

A„ = «23 • (45) 

\ / 

In essence we are proceeding as if Eq. (2) pertains to a system of two qutrits in which, however, the only wave 
functions ^„ of interest have zero projection on the eigenfunction W3 of the first qutrit. Eq. (19) now remains valid 
for this qubit-qutrit A^, but Eq. (20) must be replaced by 
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(46) 







Correspondingly, in the qubit-qutrit case the equation |det(^a)| = |det(A„)| yields merely = 0, hardly a useful 
analog of Eq. (21). 

There is a useful analog of Eq. (21), however, which is obtained as follows. From Eq. (45) 



6?i 6?2 
AaAi = I 




where 



*iil 



+ \a: 



121 + Pisl 



Prom Eq. (19), furthermore. 



7 a _ . Q* 
'^21 — "12 



22 ~ l'^2ll + 1*^221 + 1^231 ■ 



(47) 

(48) 
(49) 
(50) 
(51) 

(52) 



which expresses the singular value decomposition theorem result [6] that the diagonal elements of are the eigen- 
values of the Hermitian matrix A^A^. The corresponding eigenvalue equation implied by Eq. (52) is 



Ao^Ai - AI = C/a A2[/t - AI =C/a(A2 - \l)ul 

Eq. (53) yields 

Aiit{Ao,Al - AI) = det(A2 - AI). 
Equating the coefficients of A on the two sides of Eq. (53) we find 

^>?l^>?2-&?2&?l = (c?)'(c?)^ 



(53) 
(54) 
(55) 



a result which also follows directly from the theory [8] of the characteristic polynomial associated with A^A]^. With 
the aid of Eqs. (48)-(51) and some algebraic manipulations, Eq. (55) can be rewritten in the form 



|a?ia^2 - a?2a^i|' + K^a^g - a^ga^il' + ja^^a^g - a^ga^al^ = (c?)2(c^)2. 
Eq. (56) is the desired qubit-qutrit analog of Eq. (21). Next introduce the three matrices 



(56) 



S,, — 














1 













-1 































1 











1 

















Vo 














0/ 


/° 














1\ 





























-1 














-1 





























Vi 














0/ 



(57) 



(58) 
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(59) 



/O 0\ 

1 

-10 



-1 

\010 00 0/ 

which can be thought of as qubit-qutrit analogs of Eq. (16). In terms of these matrices Eq. (56) becomes 

1 /2 

*aS,*„|%|*„S^*„|%|*„S,*„|'^ =2c?c^, (60) 

where now is the column matrix whose six components ^'ma, rn = 1 to 6, are respectively the coefficients 
a"j^, a"2, ai3, 021' '^22' ^^23 appearing in Eqs. (2) and (45). Evidently the left side of Eq. (60), like the two qubit 
concurrence of \1/q, defined in Eq. (15), precisely equals the parameter defined in Eq. (7). Thus the left side of 
Eq. (60) can be regarded as the qubit-qutrit analog of the two qubit concurrence, and legitimately can be denoted 
by Con{'i'a). The corresponding analog of Eq. (22), expressing the minimum average concurrence any decomposition 
of p can attain, is 



JV 



N 



C(p) = min ^ paCon(*c«) = min ^ pa 



a=l 



+ 



1/2 



(61) 



Especially in view of the square roots therein, the right side of Eq. (61) is sufficiently more complicated than the right 
side of Eq. (22) that this writer has been unable to generalize Wootters' derivation to obtain a closed formula for the 
EOF of an arbitrary qubit-qutrit mixed state. A lower bound on this EOF can be derived, however, via consideration 
of the following exercise: 

We seek the minimum value fmin of 



JV 



/2 



(62) 



subject to the constraints 



N 



a=l 



N 



a=l 



N 
a=l 



z, 



(63) 



(64) 



(65) 



wherein all the Xa,ya, Za arc real variables > 0, and X, Y, Z are given fixed quantities. Then routine employment of 
Lagrange multipliers to take account of the constraints finds that the right side of Eq. (62) has the single extremum 
(which can be shown to be a minimum) 



(66) 



Eq. (66) implies 



N 



/2 



(67) 



whenever Xa,ya, Za are constrained by Eqs. (63) - (65). 
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We now observe that Eqs. (23) and (24) remain valid in the qubit-qutrit system, except that the sum over s in 
those equations now runs from 1 to 6, not merely 1 to 4; correspondingly W now is a 6x6 matrix. Similarly Eqs. (25) 
- (27) also remain fully valid, except that R and $ now are 6x matrices, remembering of course that = ^/pZ^o:■ 
In terms of the Eq. (61) takes the rather more convenient form 



TV 



C(p) = min^ I $aSa;«'a 



a=l 



1/2 



(68) 



Any given decomposition of p, i.e., any given set in the qubit-qutrit analog of Eq. (25), determines the values of 
the quantities Xd, Yd, Zn defined by 



N 



a=l 



(69) 



N 



(70) 



N 



a=l 



(71) 



$ S $ 



Identifying 
evident that Eq. (67) implies 



$ S $ 



$ s $ 



in Eqs. (69) - (71) with Xce,ya,Za respectively in Eqs. (63) - (65), it is 



N 



+ 



+ 



1/2 



where XdjYd, Zd are defined by Eqs. (69) - (71). Eqs. (68) and (72) in turn imply 

min(X|, + yJ + Z|,)V2<c(p). 
Next define for this qubit-qutrit case, in analogy with the first equality in the two qubit Eq. (28), 

N 



Cx{p) = minXD = min^ l>„Sa;$o 



a=l 



(72) 



(73) 



(74) 



N 



Cy{p) = mmYo = min^ ^'aS^^^c 



(75) 



N 



Cz{p) = minZi) = min ^ l>«S^$c 



Then Eq. (73) yields 

(C2 + c2+c2)V2<c(p). 
Furthermore, just as in Eq. (28), we now have 

N N 



C^(p) = min^ (*S^*)aa = min^ (RWS^WR)„c« 



(76) 



(77) 



(78) 
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minimized over all possible R satisfying Eq. (26). Correspondingly Eqs. (29) - (34) arc pertinent to Eq. (78) 
provided: the subscript x now is attached to the quantities C/w, A^k, A and Q; Qx now is a 6xN matrix; the sums 
over m in Eqs. (33) and (34) run from 1 to 6; and the six eigenvalues 

Axl > Ax2 > Ax3 > Ax4 > Ax5 > Ax6 > (79) 

of Awx are the square roots of the eigenvalues of p^^'^SxP*SxP^^^- 

The matrix Sx defined by Eq. (57) has two rows and two columns that arc identically zero. It readily can be seen 
that the corresponding rows and columns of the matrix Sa;/3*Sa; are identically zero as well. Thus the 6x6 Hermitian 
matrix SxP*Sx has at least two zero eigenvalues, i.e., has [11] a rank no greater than 4. It then follows [12], as can be 
verified by direct multiplication, that the Hermitian matrix p^/'^SxP*SxP^/^ also has rank no greater than 4, i.e., that 
at least A^s and Xxe arc zero in Eq. (79). Consequently, except that subscripts a; must be appropriately attached, Eq. 
(35) remains valid as written; in particular Eq. (35) gives a lower bound on Cx{p) even though it does not include 
terms proportional to Xx5 and A^e- 

Therefore, now replacing the two qubit 4x4 matrix of Eq. (37) with the unitary 6x6 matrix 



Q3> = 2 



/I 

i 

i 
i 

\0 



1 

i 

—i 
—i 







1 

—i 

i 

—i 





1 

—i 
—i 

i 












2/ 



(80) 



and obtaining the 6x6 Q3< via the corresponding replacement in Eq. (42), we conclude that 

Cx{p) = max(0,Axi - Xx2 - Ks - Xxa) 



(81) 



is the the qubit-qutrit generalization of Eq. (43) for the quantity defined by Eq. (74). Noting that each of the matrices 
Sy and S2 defined by Eqs. (58) and (59) also has two rows and two columns that are identically zero, we evidently 
can further conclude that the quantities defined by the corresponding Eqs. (75) and (76) are given by 

Cy{p) = max(0,Ayl - Xy2 - XyS - Xy4), (82) 



Cz{p) = max(0,A2i - A^2 - A23 - X^), (83) 

wherein the A^ and A^, each ordered as in Eq. (36), arc the square roots of the four largest eigenvalues (some of which 
may be zero) of the respective matrices p^/^Syp*Syp^/^ and p^^^SzP*Szp^^'^- 
Define 

Cs,{p) = {Cl + Cl + Cl)'/\ (84) 

where the quantities on the right side of Eq. (84) have the values given by Eqs. (81) - (83). Then Csbip) is a lower 
bound on the qubit-qutrit C(p) defined by Eq. (61). Correspondingly, remembering Eq. (14), the desired lower 
bound on the qubit-qutrit EOF is e[Czb{p)]- It can be seen, much as discussed for the two qubit case at the end of 
Section II. A. above, that the lower bounds given by Eqs. (81)-(83): (i) arc yielded by large ranges of 6xiV matrices 
Qx^Qj/jQz) with iV arbitrarily large, and (ii) remain valid even when the qubit-qutrit p has a rank r > 2 but less than 
the value 6 implicitly assumed in the foregoing derivations of those equations. 



III. THE QUBIT-QUDIT SYSTEM. 

Wo now turn to the derivation of the lower bound on the qubit-qudit EOF. This derivation is a straightforward 
extension of the qubit-qutrit derivation; there is no need to give as many details as were furnished heretofore. The 
analog of Eq. (45) is 
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(85) 
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wherein the dxd matrix Aa is made square by incorporating d-2 rows of zeros. The analog of Eq. (47) then is 

(86) 



... 

/l /it — ^21 ^22 ■■■ 

\ II I) 



wherein the only non- vanishing elements are &ii, J>i2, &21 ^^'^ ^221 given by: 



^?2 = E<'^2; (88) 



1^21 — "12 



(89) 



&?2 = EI«2ir (90) 

Thus AaAl^ has at most two non-vanishing eigenvalues (c")^, (02)^ which continue to obey Eqs. (5) and (55). Eqs. 
(87) -(90) now yield 



^ d d d d 

{cific^f = ^n^22 - ^?2^2i = 2 E E \« - = E E - ay<r • (91) 

i=l j=l j>i i=l 

That Eqs. (87)-(90) imply Eq. (91) can be demonstrated by straightforward algebraic manipulations, as well as via 
mathematical induction starting with the presumed correctness of Eq. (91) when the sums over i and j in Eqs. (87) 
- (91) run from 1 to d-1 only. 

Let be the column matrix whose components ^ma, m = 1 to 2d, are the coefficients ag-, in the order 

a"i, a"2, aid, afii '^221 '^2d- Define the set of d(d-l)/2 symmetric 2 rfx 2d square matrices S'^, 1 < i < d — 1, j 
>i, to be the matrices whose elements S^^ all are zero except for: 

= S^,, = 1, (92) 



= = -1- (93) 



Then 



and 



*„S'^*„ = 2(a?,a?^. - a?,.a^J, (94) 



V 1/2 

Af N / d-1 . 2 



C(p) = min^ p,Con(*„) = min^ ^^|$„S^%„ | . (95) 



Further defining 



N 



a=l 



Cij{p)=inmJ2\^aS''^a , (96) 



Eq. (77) generalizes to 
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[$:X^C^^-(p)]^<C(p). (97) 

j>i i=l 

The left side of Eq. (97) is the desired lower bound Cdb{p) on the qubit-qudit C(p), using 

Cij {p) = max(0, Ai^' - - - ) , (98) 

wherein the A*-' , ordered as in Eq. (36), are the square roots of the four largest eigenvalues of the matrix 
^i/2gy^*gij^i/2 rpj^g desired lower bound on the qubit-qudit EOF is e[Cdi,{p)]. 

IV. CONCLUDING REMARKS. 

It is needful first to point out that in a number of respects our review of Wootters' two qubit derivation [2] has 
echoed the treatment given by Audenaert et al [13]; in particular these authors have made use of the matrix notation 
for decompositions of p introduced in Eqs. (24) and (25), and have noted that the minimization over all R in Eq. 
(28) can be replaced by the more convenient minimization over all Q in Eq. (33). Also the matrices S*-' appearing in 
Eqs. (92)- (96) essentially are the "indicator matrices" introduced by Audenaert et al [13], as those matrices would be 
written for the qubit-qutrit system. For the purpose of generalizing Wootters' derivation to obtain the qubit-qutrit 
EOF lower bound, however, we have found it preferable not to rely on Thompson's Theorem [14] as Audenaert et al 
[13] do, but rather to explicitly construct the two qubit matrices Q> and Q< of Eqs. (37) and (42). This procedure 
enables our relatively straightforward derivations of the lower bound Eqs. (84) and (98), on the average qubit-qutrit 
and qubit-qudit concurrences respectively. Our procedure also makes apparent the fact that (as discussed at the end 
of Section TLB.) for any qubit-qutrit p there always actually exists a wide range of matrices Q^;, i.e., a wide range of 



which actually attains the minimum possible 



decompositions of p. having an average " ^-concurrence" 

value Cx{p) given by Eq. (81), and similarly for Cy{p) and Cz{p). 

We have performed no numerical calculations designed to estimate the utility of our lower bounds. Nor have we 
attempted, by numerical calculation or otherwise, to ascertain whether these just described three different ranges of 
decompositions of p, which separately attain (p) , Cy (p) and (p) respectively, generally will overlap sufficiently 
to ensure the attainability of the bound given by Eq. (84) . In other words we have not tried to answer the question: 
Given some arbitrary qubit-qutrtit p, will it generally be possible to find a decomposition of p whose minimum average 
concurrence equals the value of C3b(p) given by Eq. (84)? If such a decomposition exists, then [recalling the two 
qubit definition of Q given immediately preceding Eq. (32)] it must be possible to find three matrices Qx,Qy,Qz such 
that 

uLQx = ul^Qy = UI^Q, = R6, (99) 

where: Qa;,Qj;,Qz are unitary 6xiV matrices (all with the same A^) which, when employed in the qubit-qutrit analogs 
of Eqs. (33)-(35), imply Eqs. (81)-(83) respectively; Uwx,Uwy,Uwz are the properly specified analogs of the Eq. 
(29) Uw when S in Eq. (29) is replaced by Sj;,Sj,,Sz respectively (because Xx5 = Xx6 = in Eq. (79), Uwx is not 
uniquely determined by the 6x6 S,,. analog of Eq. (29), and similarly for Uwy and Uwz): and the 6x7V matrix R5 is 
the R satisfying Eq. (26) which, via the qubit-qutrit analog of Eq. (27), specifies this particular decomposition. 

Evidently it must be possible to find a set of matrices Qx,Qy,Qz satisfying Eq. (99) whenever the given qubit- 
qutrit matrix p happens to be separable, i.e., whenever E(p) = 0, in which event both C{p) and C^i,{p) also must 
vanish. On the other hand, numerical calculations by Audenaert et al [13] suggest that it may not be possible to 
find such Qx,Qy,Qz for every qubit-qutrit p with Csbip) = 0. In short, whereas each of Cx{p) = 0,Cx{p) = 0, 
Cx{p) = surely is a necessary condition for separability (as Audenaert et al [13] already had concluded, though via 
a quite different formalism than employed herein) , the totality of these three conditions may not suffice to guarantee 
separability. Moreover this writer knows of no theoretical or numerical studies which might answer the question posed 
in the preceding paragraph under the circumstance that C3b(p) > 0. In his opinion, even the demonstrated inability 
to satisfy Eq. (99) for most qubit-qutrit p with C3h(/3) ~ would not imply the inability to satisfy Eq. (99) for most 
p with C:ii,{p) > 0, especially when C:ii,{p) > results from the non-vanishing of each of C,j.(/3), Cy{p) and Cz{p). 
This opinion stems from the observation that the range of matrices Q^; yielding a specified value of Cx{p) >0 is wider 
than the range of Qx yielding Cx{p) = 0, reflecting the fact that achieving Eq. (41) requires X]m=i ^mQma = ^or 
every one of the N a— indexed columns of Q, whereas Eq. (38) can be achieved with values of X^^^i '^mQma which 
are unequal for different a, 1 < a < N. 
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Our Eqs. (56) and (91), which arc the crucial starting points for our derivations of the lower bounds C^tip) and 
Cdb(p) given by Eqs. (84) and (97) respectively, could have been deduced more elegantly but less transparently from 
the equalities 



K \ ^ K K K 



1 - Tr{p\) = [TT{pA)f - Tr{p\) = E " E '^fe = 2 E E c'^?= (10°^) 



^fe=l / fe=l l>k k=l 



dA dB dA dB 

- (aija^jamndmn ~ O-ijO-mj^mnO-in) ~ lijO-mn [o-ijO-rnn ~ '^mj^iri) ~ (100b) 

i,m— lj,n— 1 ^,m— 1 



dA dB dA dA ds ds 

— 2 ^ ^ ^ ^ l^ij^mn ^in^mjl ~ 2 ^ ^ ^ ^ ^ ^ ^ ^ \^ij(^mn ^m^mjl • (lOOc) 

i.m— lj,n— 1 myi i—1 n>j j—1 

In Eqs. (100): p = ^'a^a (elements Pij,mn = O'ijO'mn)^ is a pure state density matrix in a bipartite system composed 
of arbitrary subsystems A and B (not merely a system composed of a qubit subsystem A and a qudit subsystem 
B); the wave function ^'^ has the expansion given by Eq. (2), except that i and j now run from 1 to dA and ds 
respectively, where ds are the respective dimensionalities (i.e., the maximum number of independent orthonormal 
eigenfunctions) in subsystems A, B; pA (elements {pA)im = X]j=i '^y ^mj)' reduced density matrix associated 

with subsystem A; the K real quantities c| arc the non- vanishing eigenvalues of pA, the here superfluous index a on 
the quantities afj and of Eqs. (2) and (4) has been dropped: and the elements {pA jim for a qubit-qudit system 
can be seen to be identical with the elements hm of Eqs. (48)-(51). For the qubit-qudit system, where K = 2, the 
equality between the last term in Eq. (100a) and the last term in Eq. (100c) is Eq. (91). 

Eqs. (100), which in essence have been stated by Rungta et al [15] and by Albeverio and Fei [16], are the basis 
for those authors' proposed generalizations of Wootters' [2] two qubit concurrence. In particular, except possibly 
for here inconsequential numerical factors, the left side of Eq. (100a) is the square of the "I-concurrencc" defined 
by Rungta et al [15], while the right side of Eq. (100c) is the square of the "generalized concurrence" introduced 
by Albeverio and Fei [16]; the right side of Eq. (100c) also is the "length" squared of the "concurrence vector" 
introduced by Audenaert et al [13], again except for a here inconsequential numerical factor. The numerical equality 
of these different proposed generalizations of the two qubit concurrence, remarked on by Wootters [1], is obvious from 
Eqs. (100). The present paper shows that these generalizations are appropriate and useful in qubit-qudit systems, 
where the von Neumann entropy of a pure state is a function of a single real parameter only. On the other hand, the 
presence in Eq. (100a) of terms in c^cf beyond cf c| when K >2 strongly suggests that the utility of the aforementioned 
concurrence generalizations for estimating the EOF in bipartite systems not containing a qubit will be significantly 
less than in the K = 2 qubit-qudit systems which are the subject of this paper. 

Ultimately the analytical difficulties engendered by the square roots in Eq. (61), remarked on in connection with 
that equation, stem from the decision that the right sides of Eqs. (8) and (9) arc to be regarded as fimctions of 
Cct rather than as functions of C^, i.e., as functions of the concurrence rather than the square of the concurrence. 
Focusing on C^, would have led to replacement of Eq. (61) by 

/ 2 2 2\ 

t(p) =minEl'« ( *aS^*a + ^aSy^oc +*aS^*a ), (101) 

a=l ^ ^ 

which avoids the troublesome square roots. The quantity t{p), here denoting the minimum average square concurrence 
that can be attained by any decomposition of p, is known as the "tangle" [17]. Unfortunately this writer knows of 
no way to express the tangle of an arbitrary qubit-qutrit p in the form of readily computable properties of matrices 
simply related to p, e.g., in terms of the eigenvalues of p^^'^SxP*SxP^^^ , as in Eqs. (81)-(83) yielding the lower bound 
Eq. (84). Osborne [17] has found a closed form expression for the tangle of any bipartite density matrix of rank 
r = 2, and [noting that the von Neumann entropy of Eq. (6) is a monotonically increasing concave (not convex) 
function of via Eqs. (8)-(9)], has used this expression to deduce a very close upper bound for the EOF of a rank 
2 qubit-qudit p. We have put no rank limitations on our qubit-qudit density matrices. 

Finally, it hardly is necessary to note that any of the foregoing remarks which were made solely in a qubit-qutrit 
context have obvious generalizations to the qubit-qudit system. 
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